Minigap, Parity Effect and Persistent Currents in SNS Nanorings 
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We have evaluated a proximity-induced minigap in the density of states (DOS) of SNS junctions 
and SNS nanorings at an arbitrary concentration of non-magnetic impurities. We have demonstrated 
that an isotropic energy minigap in the electron spectrum opens up already at arbitrarily weak 
disorder, while angle resolved DOS at higher energies can remain strongly anisotropic. The minigap 
value Sg can be tuned by passing a supercurrent through an SNS junction or by applying a magnetic 
flux $ to an SNS ring. A non-monotonous dependence of Eg on <& has been found at weak disorder. 
We have also studied persistent currents in isolated SNS nanorings. For odd number of electrons 
in the ring we have found a non-trivial current-phase (current-flux) relation which - at relatively 
high disorder - may lead to a 7r-junction state and spontaneous currents in the ground state of the 
system. 



I. INTRODUCTION 

In hybrid structures composed of a superconductor 
(S) and a normal metal (N) Cooper pairs can penetrate 
into the latter thereby significantly changing the prop- 
erties of the systemi. As a result of this proximity ef- 
fect the N-metal also acquires superconducting properties 
being able to carry supercurrent and, hence, exhibiting 
JosephsonSi^ and Meissneri^ effects. 

Another interesting consequence of this proximity- 
induced superconductivity is the existence of a minigap 
in the electron spectrum of a normal metal. In the 
diffusive limit this minigap Eg was found to be of or- 
der of the Thouless energy exh of this metal£i2i2ii£. In 
strictly ballistic SN and SNS systems this gap is strongly 
anisotropic^^, and it vanishes for electrons propagating 
parallel to the SN interface. For arbitrary concentration 
of non-magnetic impurities the minigap dependence on 
the electron elastic mean free path was studied in Ref.H^ 
In SNS structures the proximity-induced minigap Sg can 
be tuned by applying the phase difference (p across the 
N-metal2iH. 

Modifications of the normal metal density of states 
(DOS) due to the proximity effect can be - and were - 
studied experimentally with the aid of tunneling spec- 
troscopy methods, see, e.g., Ref. In this paper 
we point out that the proximity-induced minigap in SN 
sandwiches can be directly measured with the aid of the 
superconducting parity effec t -^^d^i^'' . Indeed, the ground 
state energies of isolated superconducting system with 
odd and even electron numbers should differ exactly by 
the value of the gap in the electron spectrum. In the case 
of superconducting grains it is just the BCS gap, while 
in SN and SNS structures this value should be set by 
the proximity-induced minigap Eg. Hence, the latter can 
be directly observed in an experimental setup similar to 
that used, e.g., in Ref. 

The possibility of tuning the minigap value Sg^tp) by 
passing the supercurrent through SNS junctions provides 
additional ways to independently test theoretical predic- 
tions. One can, for instance, consider an isolated su- 



FIG. 1: A superconducting ring with embedded normal metal 
of length d. 



perconducting ring with an embedded layer of a normal 
metal, the so-called SNS ringiS. This system is depicted 
in Fig. 1. Applying an external magnetic flux $ to such 
a system one induces persistent currents (PC) circulating 
inside the ring. Both the magnitude and the flux depen- 
dence of such currents will depend on the parity of the 
total number of electrons in the ring^^'^". The difference 
between PC values for odd (/q) and even (/g) ensembles 
is related to the minigap value £g{^). This relation ac- 
quires a particularly simple form in the limit T — > in 
which case one findsi^ 



/„($) =/e($) + 2e 



d(p 



(1) 



The last term in this equation describes the contribu- 
tion to the current from the "odd" electron occupying 
the lowest available state above the minigap £g{(p) in the 
density of states of the normal metal. 

For a broad range of system parameters the phase dif- 
ference (p across the SNS junction is linked to the exter- 
nal magnetic flux $ by means of the standard relation 
ip = 27r$/$o ('i'o is the flux quantum) which will also be 
assumed to hold throughout this paper. Eq. can be 
used, on one hand, for an independent study of the phase 
dependence of the minigap Sg{ip) and, on the other hand, 
for further investigations of the parity-affected persistent 
currents in SNS nanorings. 
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The structure of our paper is as follows. In Sec. II 
we will employ the quasiclassical formalism of Eilen- 
berger equations and evaluate the phase-dependent mini- 
gap £g{^p) in SNS systems at an arbitrary concentration 
of non-magnetic impurities. We will demonstrate that an 
isotropic minigap opens up in the normal metal already 
in the limit of very weak disorder. At sufficiently large 
values of the electron elastic mean free paths we recover 
a non-monotonous dependence of the minigap on ip. In 
Sec. Ill we will analyze an interplay between parity ef- 
fect and persistent currents in SNS nanorings with an 
arbitrary impurity concentration. We will show that in 
the case of odd total number of electrons in the ring the 
dependence of PC on the applied magnetic flux acquires 
non-trivial features which might lead to a ir-junction be- 
havior and spontaneous currents in the ground state of 
the ring. A brief discussion of our key observations is 
presented in Sec. IV. 



II. PHASE-DEPENDENT MINIGAP AT 
ARBITRARY IMPURITY CONCENTRATIONS 

Let us analyze thermodynamic properties of SNS rings 
presented in Fig.^ We will stick to the simplest case of a 
superconductor with singlet isotropic pairing and assume 
that both NS interfaces are fully transparent. The length 
of the normal metal layer and its cross section are denoted 
respectively as d and A. 

Our analysis is based on the Eilenberger equationsiSiSi 
for the energy-integrated retarded 2x2 matrix Green 
functions g 



efa - A(r) - t{r,£),g{pp,r,£) = ~ivpyg{pp,r,£), 

(2) 

f{pF,r,e) = l, (3) 



where [a, b] = ab — ba, e is the quasiparticle energy, pp = 
mvp is the electron Fermi momentum vector and fa is the 
Pauli matrix. The matrices g and A have the standard 
form 
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where A is the BCS superconducting order parameter 
chosen to be spatially constant in the superconductor and 
equal to zero in the normal metal. Electron scattering 
on non-magnetic impurities is accounted for by the self- 
energy term S. Within the Born approximation one has 



t{r,e) = {g{pp,r,e)), 



(5) 



where £ is the electron elastic mean free path and angular 
brackets (...) denote averaging over the Fermi momentum 
directions. 
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FIG. 2: Angle averaged DOS in the middle of the normal 
layer. Here and below we choose d — lO^o. 



The current density and the angle resolved local den- 
sity of states are defined as 

oo 

j{r) = eNoRe J rfetanh (^^^ {vpg{pp,r, e)) , (6) 

— OO 

iy{pF,r,e) :^Reg{pp,r,e), (7) 

where DOS i'{pp,r,e) is normalized to its normal state 
value at the Fermi energy A^o — mpp /2tt^. Let us 
also note that the above quasiclassical formalism and, 
in particular, the description of electron scattering in 
terms of the self-energy ^ are applicable provided the 
number of conducting channels in the system is large 
Nch = App/in :g> 1. This inequality will be assumed 
to hold throughout the whole calculation^^. In addi- 
tion, we will assume that the length of the N-metal 
strongly exceeds the superconducting coherence length 
d » ^0 ~ VP /A. 

Our first goal is to analyze the behavior of DOS 
^{pp^r^e) in the normal metal as a function of the phase 
difference tp = 27r$/$o for arbitrary values of elastic 
mean free path t. To begin with, let us recall that in 
a strictly ballistic limit € ^ oo the spectrum of an SNS 
junction below the BCS energy gap is formed by discrete 
Andreev levels with energies En- In the limit En <IC A 
the values En readSiii 



En = ^[^^(2"- + 1) + '/ssgnwj;] 



(8) 



where Vx is the x-component of the Fermi velocity vector 
vp. This expression demonstrates that the proximity 
induced gap in the spectrum of a ballistic SNS system 
is strongly anisotropic. It depends on the direction of 
the Fermi velocity and vanishes for electrons propagating 
parallel to NS interfaces. In other words, the trajectories 
for such electrons do not cross NS interfaces and, hence, 
these electrons do not suffer Andreev reflection and do 
not "feel" any proximity effect. 

The situation changes significantly in the presence of 
already very weak disorder. In this case all electrons 
in the N-metal get scattered by non-magnetic impuri- 
ties and, sooner or later, their trajectories hit one of the 
NS interfaces and electrons get Andreev-reflected. As a 
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FIG. 3: The angle resolved DOS in the middle of the normal 
layer (thin lines) together with the angle averaged DOS (bold 
solid line). The angle between quasiparticle momentum and 
the x-direction (normal to NS interfaces) is denoted as 9. 



result, an isotropic proximity-induced minigap £g{^) de- 
velops in the density of states of the normal metal. The 
value of this minigap depends both on the electron mean 
free path £ and the phase difference Lp but it remains in- 
dependent of the momentum direction for all values of ^ 
and Lp. 

This property is illustrated in Figs. El and |31 Fig. 
121 shows DOS inside the normal metal with disorder to- 
gether with that evaluated in the ballistic limit ^ oo. 
An important difference between the two curves is ob- 
served at low energies: The proximity-induced minigap 
is strictly zero in the ballistic limit while it opens up 
in the presence of disorder. At higher energies electron 
scattering on disorder broadens and eventually washes 
out peaks corresponding to Andreev levels © . 

Fig. displays a typical energy dependence of the an- 
gle resolved density of states v{pp,r, e) inside the normal 
layer of an SNS structure in the presence of disorder. We 
observe that v{pF,r,e) vanishes for all e < £9(1^) in- 
dependently of the direction of the Fermi momentum. 
Thus, already very weak disordes^i makes the minigap 
isotropic for all values of Lp. In addition, we have verified 
that the minigap value does not depend on the coordinate 
inside the normal metal. 

At the same time at energies just above the minigap 
s > Eg the anisotropy in the density of states is clearly ob- 
servable for sufficiently large values of the electron mean 
free path i, see Fig. This anisotropy decreases with 
decreasing mean free path £, and in the dirty limit i ^0 
the density of states v[pf,r, e) becomes almost isotropic. 
Anisotropic behavior of the density of states in the ultra- 
clean limit £ ^ d can easily be understood if one ob- 
serves that the minigap value in the latter limit is of 
order Eg ~ u_f/^ ^ vp/d. Hence, for quasiparticles with 
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FIG. 4: The minigap £5(0) and the zero temperature Joseph- 
son critical current of an SNS junction as a function of d/l. 
For sufficiently short values of I both quantities approach the 
corresponding results derived in the diffusive limit from the 
Usadel equation analysis. 



energies just above the minigap vp/d » e > vp/£ one 
can estimate Vx ~ vpdji <C w_f, i.e. these states are 
formed predominantly by grazing electrons. 

Let us also note that in the normal metal at energies 
just above the minigap e > Eg the angle averaged DOS 
v{r,E) = {i^{pp,r,E)) depends on energy as 



j/(r,e) cx (e - Eg) 



1/2 



(9) 



This square-root dependence appears to be universal for 
all mean free path values £. In order to verify this state- 
ment we have plotted In i/(e) versus logarithm of the pa- 
rameter (e — Eg) /Eg for different values of £ (not shown). 
For all impurity concentrations at energies above the 
minigap our numerical data points collapse on straight 
lines corresponding to the same slope 1/2. For very clean 
systems the dependence © is observed within narrower 
energy interval above the gap, and in the limit £ 00 
this interval shrinks to zero faster than the minigap itself. 
On the contrary, for dirtier systems Eq. Q applies for 
broader energy intervals and match with the correspond- 
ing dependence^ established from Usadel equations in 
the diffusive limit. 

The dependence of the minigap Eg at ip — on the 
mean free path £ is depicted in Fig. 01 We have evaluated 
the minigap at d = lO^o and d = lO'^^o ^^nd found a very 
good agreement with earlier numerical results of Pilgram 
et alm^ obtained for NS structures with the normal layer 
thickness equal to lO^^o- Our numerical results for Eg = 
Eg (0) are rather well approximated by a simple analytical 
formula 



Vp 



ay 



d hy"^ 



cy - 



1' 



where y = dl£ and a « 0.47, b k, 0.45, c « 
decreasing £ the minigap first grows as Eg 

Q.nvp/d at £ 



(10) 

1.45. With 
a QA7vp/£, 
a 0.67d and 



reaches its maximum e™'^'^ 
then decays for small £ as Eg ~ 3.12eTh = l-04:Vp£/d^ 
approaching the results of Belzig et aZAiS derived from 
the Usadel equations analysis. Here and below ETh — 
vp£/3d'^ is the Thouless energy of a normal metal. 

Phase dependence of the minigap is shown in Fig.[Slfor 
the mean free path values ranging between quasi-ballistic 
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FIG. 5: Phase dependence of the minigap eg{(f) for an SNS 
structure at different values of the electron mean free path I. 



[l ^ d) and diffusive (I ^ c?) regimes. A somewhat 
unexpected non-monotonous dependence of the minigap 
value on the phase difference is observed in the ultra- 
clean limit £ ^ d. We believe that this effect is caused 
by a non-trivial redistribution in the angle resolved DOS 
in the presence of the phase twist cp across the N-layer. 
With decreasing the mean free path the maximum of the 
function eg{ip) moves away from the point ip = n to- 
wards smaller values of the phase difference. Eventually 
the maximum of Sg^tp) reaches the point ip = and this 
dependence becomes monotonous at all smaller values of 
£. In the diffusive limit I <^ d our results approach those 
obtained by Zhou et al}^ from Usadel equations, cf. Fig. 
2 in Ref. 13. 

In the limit £ ^ d the minigap dependence on p is 
rather well approximated by a simple parabola 



(11) 



Note that although this dependence is in agreement with 
our numerical results as well as with those displayed in 
Fig. 2 of Ref. 13, it does not agree with Eq. (5) from 
that paper at small phases <C tt. 



III. PERSISTENT CURRENTS AND PARITY 
EFFECT 

We now go over to the calculation of persistent cur- 
rents in isolated SNS nanorings. As the total number of 
electrons in the ring is fixed to be either even or odd 
it is necessary to account for the parity effect. This 
task will be accomplished within the parity projection 
technique2i2SiS& recently adaptedi^iSS to the calculation 
of PC in superconducting nanorings. 

Let define and evaluate the even/odd parity projected 



thermodynamic potentials 

,'1 ± e"^"! 

r^e/o = rif — Tin 



/ 



A = {nh-nf)/T, (12) 



where fiy is the standard grand canonical thermody- 
namic potential and Qb is obtained from ri/ by expressing 
the latter as a sum over the Fermi Matsubara frequencies 
ujf = 7rr(2n -|- 1) with subsequent substitution of ujf by 
the Bose frequencies ujb — 2'KTn. For further details we 
refer the reader to the paper—. 

With the aid of the Eilenberger quasiclassical 
formalismSi it is easy to cast the grand canonical ther- 
modynamic potential of the system fi/ to the following 
form 

oo 

rif^n-moT j dr j de\n[2cos\i{e/2T)]v{r,e), (13) 



where = / dr|A^(r)|/A + X]fc ^fe^ ^ is the BCS coupling 
constant and S^k is the single particle energy. Similarly, 
for rj;, we find 



nb = n-ANaT j dr j rfe In [2 sinh(e/2r)] zy(r, e). (14) 



With the aid of the above expressions it is now possible 
to directly evaluate PC circulating in a superconducting 
ring both for even and odd electron ensembles. Taking 
the derivative oi^j /b with respect to the phase difference 



T T + ^"-^Z 

^ exp(A)±l- 



(15) 



where the upper (lower) sign corresponds to the even 
(odd) ensemble and 



/o 



dflf/b 



(16) 



It is well known that the parity effect is mostly pro- 
nounced in the low temperature limit in which case the 
parameter A H12I) is much larger than oneS^ 1. Eval- 
uating the integrals in Eqs. H13I) . 114(1 at temperatures 
T ^ |eg('y3)| we obtain 

AiT)^AMiT)+AsiT), (17) 

where 

An{T) ^ NaVi,T^/\-^'^{^) eM-sM/T). (18) 

As{T) - A^ol^sT^/^Ai/^ exp(-A/T) (19) 

are the contributions respectively from the normal and 
superconducting parts of the ring with the corresponding 
volumes Vn and Vs- The term Aj^ was evaluated making 
use of the dependence Q for DOS at energies e just above 
the minigap. 
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FIG. 6: Phase dependence for the zero temperature Josephson 
current at different values of the mean free path I. 



Since throughout this paper we always assume the 
minigap to be smaller than the superconducting gap, 
eg{(p) < A, at sufficiently low T the term An domi- 
nates over As provided the volume ratio Vn /Vs is not 
too small. Assuming that the contribution As can be 
neglected, from the condition^^ A{T*) 1 we arrive at 
the estimate for the crossover temperature T*. above 
which the difference between thermodynamic potentials 
for even and odd ensembles is negligible and the parity 
effect is smeared out by thermal fluctuations. Within the 
logarithmic accuracy we flnd 



T* ^egiip)/HNaVNeg{ip)). 



(20) 



In the whole temperature range x < T* where the parity 
effect remains pronounced the temperature effect on the 
Josephson current across SNS junction turns out to be 
negligible for all values of the mean free paths I. This is 
guaranteed by the inequality T* <^ min(£Th, vp/d) which 
is always satisfled for generic systems. Hence, in order to 
study PC in the presence of the parity effect it is suffi- 
cient to restrict our calculation to the zero temperature 
limit T = 0. In this case from the above analysis one 
immediately arrives at Eq. which establishes a di- 
rect relation between PC values in the even and the odd 
superconducting ensembles. 

Let us first evaluate PC for the even ensemble /g. At 
T = this current identically coincides^^'^^ with one cal- 
culated for the grand canonical ensemble. The latter is 
easily evaluated with the aid of the quasiclassical Eilen- 
berger equations H2I6|I . The results for the current-phase 
relation /e(v) displayed in Fig.|Slfor various impurity 
concentrations. The dependence of the critical current Ic 
on the electron mean free path £ is presented in Fig. 0] 

PC in the odd ensemble /q at T = can now be triv- 
ially evaluated by making use of Eq. ^ and combining 
our results for Ie{^) with those for the minigap Sg{ip) 
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FIG. 7: Phase dependence of the Josephson current at T = 
for the odd and even number of electrons in the ring. 



derived in the previous section. The typical dependence 
/o(iy9) is displayed in Fig. [7| We observe that at suffi- 
ciently large values of 93 < tt the absolute value of the 
odd electron contribution to PC ledegjdip exceeds the 
term /e(v) and the total current /q changes the sign. 
This non-trivial parity-affected current-phase relation is 
specific for SNS rings with disorder and it substantially 
differs from the current-phase relations derived earlier for 
SNS rings with ballistic^^ and resonantSS transmissions. 

At the same time, as in the previous case^i^iSli, in the 
odd ensemble there emerges a possibility for a 7r-junction 
state as well as for spontaneous currents in the ground 
state of the system without any externally applied mag- 
netic flux. In order to analyze the situation it is sufficient 
to evaluate the ground state energy of the SNS junction 
by integrating Eq. with respect to the phase ip. One 
finds 





(21) 

where E^ig{tf) are the ground state energies of SNS 
junction for even and odd number of electrons in the 
ring. While the energy Ee{ip) is always non- negative and 
reaches its minimum at 1^9 = 0, in the odd case the ground 
state energy Eo{^p) can become negative reaching its ab- 
solute minimum at (p = tt. This physical situation of a 
TT-junction is illustrated in Fig. [S] 

It is easy to find out under which conditions the tt- 
junction state becomes possible. For that purpose it is 
sufficient to observe that for any impurity concentration 
£'e(7r) = a/c/e, where Ic is the grand canonical critical 
current at T = and a is a number of order one which 
depends on the particular form of the current-phase re- 
lation. The TT-junction condition Eo{'n) < is equivalent 
to the inequality 



eg(0) > alc/e. 



(22) 



From Fig. ^ it is obvious that in the many channel 
limit the inequality (|22l) cannot be satisfied for suffi- 
ciently large £, in which case Ic is large and, on the 
contrary, the minigap £^(0) is smallSi. On the other 
hand, for sufficiently short values of the mean free path 
Ic oc P decays faster with decreasing t as compared 
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FIG. 8: Josephson energy E{ip) of an SNS ring as a function 
of the phase difference ip for the even and odd ensembles. The 
soUd curve corresponds to a vr-junction state. 



to the minigap £g(0) oc £, and the Tr-junction state be- 
comes possible. In particular, in the diffusive limit one 
findgio Ic ^ 10.82£Th/ei?w = l.bSeNchVpi^ /cP and 
a ~ 1.05, where i?Ar is the Drude resistance of a nor- 
mal metal. Combining these results with the expression 
for the minigapSii^ £g{0) — 3.12eTh, from the condition 
(|22|l we observe that in the odd case a 7r-junction state is 
realized provided the number of conducting channels in 
the junction Nch is smaller than 



Nch < OMd/e. 



(23) 



This condition is not very restrictive and it can certainly 
be achieved in various experiments. For sufficiently dirty 
junctions it allows for a formation of a 7r-junction state 
even in the many channel limit. The condition H23|) can 
also be rewritten as 

gN < 1.73, 

where gN — Rq/Rjq — 8Nchl/3d is the dimensionless 
conductance of a normal layer and Rq — h/e^ k, 25.8 kil 
is the quantum resistance unit. 

The condition for the presence of spontaneous currents 
in the ground state of SNS rings with an odd number of 
electrons is established analogously, one should only take 
into account an additional energy of the magnetic field 
produced by PC circulating inside the ring. The ground 
state with spontaneous currents is possible provided the 
total energy of the ring -EtotlTr) becomes negative, i.e. 



StotW = 1.8eTh[5iv-1.73] 



(^o/2)^ 
2C 



< 0, (24) 



where C is the ring inductance. This condition is more 
stringent than that for the 7r-junction state, but can also 
be satisfied provided the ring inductance C exceeds a 
certain threshold value which, strictly speaking, depends 
on gpf and can roughly be estimated as ~ 0.1$Q/eTh- 



IV. DISCUSSION 

For clarity, let us briefly summarize our key observa- 
tions. Our analysis was focused on the two main issues. 



the proximity-induced minigap in NS and SNS structures 
at an arbitrary concentration of non-magnetic impurities 
and the effect of parity number on persistent currents 
in SNS nanorings. We have demonstrated that already 
weak disorder qualitatively modifies the density of states 
in the normal metal of NS and SNS proximity structures. 
A striking observation at this point is that the proximity- 
induced minigap Eg turns out to be isotropic for an ar- 
bitrarily weak disorder—, even though DOS at energies 
above the gap £ > £g may remain highly anisotropic, cf., 
e.g., Fig. O Another interesting observation is the possi- 
bility of a non-monotonous dependence of the minigap Eg 
on the applied magnetic flux $ (or phase i^) in SNS rings, 
cf. Fig. El Such a non-monotonous dependence can be 
realized only in the limit of weak disorder, whereas for 
not too long elastic mean free paths a monotonous de- 
crease of the minigap from its maximum value at $ = 
to zero at $ = <&o/2 was found. 

Piercing an SNS ring by an external magnetic flux one 
induces circulating persistent currents in such a ring. In 
an isolated ring both the amplitude and flux dependence 
of such PC may strongly depend on the electron parity 
number. Provided the number of electrons in the ring is 
odd, one electron remains unpaired down to T = and 
occupies the lowest available state above the proximity- 
induced minigap Eg{(p). This electron produces a coun- 
tercurrent which - in the limit of relatively short electron 
mean free paths - may signiflcantly modify the current- 
phase relation and yield a 7r-junction behavior and spon- 
taneous currents in the ground state of an SNS ring flow- 
ing in the absence of any external magnetic flux. Note, 
that although these observations look qualitatively simi- 
lar to earlier results derived for SNS rings with ballistic^S 
and resonant^'' transmission, there also exist important 
differences. In particular, the current-phase relation is 
entirely different in the diffusive limit considered here. 
Also the restriction on the number of conducting chan- 
nels Nch in the normal metal \2'S\i is less stringent that 
that formulated in Refs. This feature of diffusive 

SNS rings is rather advantageous for possible experimen- 
tal observation of the effects discussed here. 

Our results demonstrate that superconducting parity 
effect can be used in order to directly measure both the 
magnitude and the flux dependence of the minigap in 
SNS nanorings. This can be done, e.g., with the aid of a 
setup similar to one used in Ref. [T3. Such type of mea- 
surements is not restricted by the number of channels in 
the normal part of the ring and may serve as an alterna- 
tive to tunneling spectroscopy of the minigap. Supercon- 
ducting rings with embedded carbon nanotubes might be 
promising candidates in order to experimentally investi- 
gate parity-affected persistent currents in SNS nanorings. 
Several gj-oups^^i^SiSiiS have recently reported observa- 
tions of dc Josephson current in superconducting junc- 
tions with carbon nanotubes. Therefore, it appears feasi- 
ble to fabricate and experimentally investigate SNS rings 
with carbon nanotubes which should exhibit properties 
predicted and analyzed in our paper. 
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